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REFIECTION OF A WEAK SHOCK WAVE FROM A BOUNDARY 
LAYER ALONG A FLAT PLATE 

I - INTERACTION OF WEAK SHOCK WAVES 
WITH LAMINAR AND TURBULENT BOUNDARY LAYERS 
ANALYZED BY MOMENTUM- INTEGRAL METHOD 
By Alfred Ritter and Yung -Huai Kuo 

SUMMARY 


The present paper is concerned with the phenomena encountered when 
a plane oblique shock wave is incident upon the boundary layer of a flat 
plate. In an effort to simplify the problem, the flow field was divided 
into a viscous layer near the wall and a supersonic potential outer flow. 
The pressure disturbances due to the incident wave would be propagated 
upstream and downstream in the subsonic portion of the boundary layer, 
thus giving rise to perturbations of the boundary layer. By restricting 
the study to infinitesimal incident compression waves, only small per- 
turbations were encountered and hence the ordinary linearized theory 
could be applied to the outer flow. In the laminar case, the boundary- 
layer treatment was based upon a momentum- integral equation previously 
derived by Howarth. The two flows must be compatible; hence, the deflec- 
tion of the stre amlin es near the boundary layer was expressed in terms 
of the vertical velocity component along the edge of the boundary layer 
and this relation was used as a boundary condition for the outer flow. 

The boundary condition determined the form of solution upstream and down- 
stream of the point of incidence. Determination of the constants of 
integration was accomplished by a consideration of conditions at infinity 
and a matching of the two flows at the point of incidence. With the 
outer flow thus determined, boundary-layer growth and pressure distribu- 
tion were computed and results for the 1 ami nar case were obtained as 
follows : 

(a) The pressure disturbance along the wall decreased exponentially 
from a definite value at the point of incidence to zero far upstream of 
the point of incidence. Downstream of the point of incidence, the pres- 
sure rose to a maximum value and then dropped off to the value corre- 
sponding to regular reflection. 
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(b) The disturbances produced by the interaction decayed exponen- 
tially upstream; for a free -stream Mach number of approximately 2 and a 
Reynolds number of approximately 1500 in the undisturbed boundary-layer 
displacement thickness the upstream influence was of the order of 30 
boundary-layer displacement thicknesses. 

(c) The "self-induced" pressure gradient along the wall was such 
that the boundary layer might separate ahead of the point of incidence. 

If separation occurred, the separation point moved upstream as the shock 
strength was increased. With increasing Reynolds number, the separation 
point also moved upstream, whereas for increasing Mach number, the sepa- 
ration point moved downstream. 

In the turbulent case the upstream influence was quite small and the 
incident wave must be reflected as a shock wave. 


IMRODUCTION 


It has been found that if the free-stream subsonic Mach number 
becomes high enough so that local supersonic zones are formed on an air- 
foil, sharp changes in the airfoil characteristics occur which cannot be 
explained by classical aerodynamics. The resulting loss of lift is 
accompanied by a large increase of drag in consequence of the appearance 
of shock waves on the surface of the airfoil. However, close study 
reveals that the drag increase is too large to be accounted for by the 
shock loss and change in skin friction. It must be caused by the sudden 
change of the flow pattern. This seems to indicate that the shock wave, 
when formed over the airfoil surface, modifies the character of the 
boundary layer in such a way as to create a wider wake. 

At Guidonia, Ferri, by examining the measured pressure distributions 
over airfoils at supersonic speeds, found that for the forward portions 
of the airfoil the experimental pressure distribution agrees quite well 
with that calculated from potential-flow theory (reference l) . As the 
trailing edge is approached, however, the experimental pressures along 
the upper surface, for positive angles of attack, become considerably 
higher than the calculated ones. (At negative angles of attack, this 
behavior occurs along the\Lower surface.) This discrepancy is under- 
standable in the light of the fact that the boundary layer on the rear 
portion of the airfoil is extremely sensitive to pressure disturbances . 
Since the flow far away from the airfoil has to return to its original 
pressure and direction, a shock wave must emanate at or near the trailing 
edge. There is a sharp pressure rise across the shock and pressure dis- 
turbances are transmitted upstream in the subsonic portion of the boundary 
layer causing a thickening of the layer near the trailing edge. This 
thickening, in turn, generates compression waves which travel downstream 
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and interact with the shock wave. The net result is that there is, 
starting from some point forward of the trailing edge, a gradual compres- 
sion to the main-stream pressure. The pressure distribution is hence 
altered in such a manner as to give values of lift and pressure drag 
smaller t h a n those calculated from theory. The measured moment coeffi- 
cient will also differ from the theoretical values ■ 

The preceding discussion shows that, in the transonic and supersonic 
regimes, the Mach number alone is insufficient to deter min e the flow 
characteristic sj the Reynolds number can also be important. This was 
first demonstrated by Ackeret, Feldmann, and Rott in .experiments which 
established the close relationship between the shock-wave pattern and 
the Reynolds number (reference 2). These prove beyond doubt that the 
flow far away from the wall depends intimately upon the character of the 
boundary layer, that is, the Reynolds number. Thus, under such circum- 
stances, the concept of the boundary layer requires modification. 

As a first step, a simplified problem of a plane shock wave incident 
upon the laminar boundary layer over a flat surface will be dealt with. 
For a number of years. Dr. H. W. Liepmann and his associates at the 
California Institute of Technology have been conducting experimental 
studies of shock-wave boundary-layer interaction. Recently they investi- 
gated the problem of the reflection of shock waves from boundary layers 
(reference 3)- Based on the qualitative experimental data, some impor- 
tant conclusions were drawn regarding the reflection of an incident shock 
wave from a boundary layer: 

(a) _ The type of boundary layer, whether laminar or turbulent, 
markedly affects the interaction. With a turbulent boundary layer, the 
reflection is practically the same as that for a regular reflection in 
nonviscous flow theory. For a laminar boundary layer, however, there is 
a large interaction zone near the point of incidence that is quite dif- 
ferent from the regular invlscid reflection. The incident wave is appar- 
ently reflected as from a constant-pressure surface in the form of a 
Prandtl -Meyer fan. In returning to its final direction parallel to the 
wall, the deflected flow, behind the expansion, recompresses to the pres- 
sure appropriate to that behind a regular reflection. 

(b) For M = 1.4 and Re x = ^ = 0.9 X 10^, the influence of the 

Vqo 

incident wave extends upstream in a laminar layer a distance of the order 
of 50 boundary-layer displacement thicknesses. On the other hand, the 
upstream influence is practically negligible for the turbulent boundary 
layer. 

(c) It appears that, except for very weak incident waves, the laminar 
layer almost always separates somewhere in front of the point of inci- 
dence, whereas no separation was found for the turbulent layer. 
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Quite recently, Barry, Shapiro, and Neumann (reference 4) made a 
similar study and their results are in substantial agreement with the 
GAICIT results. An advantage of their experiments, however, is that, 
although the tests were carried out at constant Mach number, the Reynolds 
number and shock strength were varied so that more qualitative as well as 
q uant itative data are available to serve as a guide to future theoretical 
work. 


Theoretical solutions of this. problem have been almost as meager as 
the experimental results. A first attempt was made by Howarth (refer- 
ence 5), who considered the case of a wave incident upon the interface 
bounding two semi-infinite uniform streams, one supersonic and one sub- 
sonic. Viscosity and heat conduction were neglected and the equations 
were lin earized - He then demonstrated the upstream propagation of dis- 
turbances in the subsonic portion of the flow field and showed that an 
incident compression wave is reflected as compression upstream of the 
point of incidence -while being reflected as expansion downstream of the 
point of incidence. However, these results can only be regarded as quali- 
tative because of the obvious shortcomings of the model. In an effort 
to make the Howarth model more realistic, Tsien and Finston simulated 
the boundary layer with a uniform subsonic stream bounded by a wall on 
one side and a semi-infinite, uniform, supersonic stream on the other 
(reference 6). As before, viscosity and heat conduction were neglected 
nnd only small disturbances were considered. For the case of a compres- 
sion wave incident upon the interface separating the two streams, it was 
shown that except for the local interaction, the incident compression 
wave is regularly reflected. Locally, however, the pressure along the 
interface exhibits compression ahead of the point of incidence and expan- 
sion immediately behind it. This local condition is qualitatively the 
same as that observed in experiments with Bhock reflection from laminar 
boundary layers} consequently there arose the speculation that perhaps 
the effects of viscosity and heat conduction were actually not too impor- 
tant in comparison with the effect of coexistence of supersonic and sub- 
sonic streams. If this were the case, it should follow that an improve - 
ment in the Tsien-Finston' model could give rise to more quantitative 
results. Now an obvious discrepancy between the Tsien-Finston model and 
the physical, situation is that the uniform subsonic stream, in attempting 
to s imula te the bcrundary layer, is incapable of satisfying the "no-slip" 
condition at the wall. Thus, the next logical step in the development 
of the theory would appear to be that of considering a main-stream velo- 
city that varies from zero at the wall to a given uniform supersonic 
velocity a short distance from the wall, and, in fact, this has been con- 
sidered by Robinson (reference j) and Lighthill (reference 8). Robinson 
assumed that the main-stream velocity varies continuously from zero at 
the wall to supersonic velocity some distance from the wall, whereas 
Lighthill assumed that the Mach number varies from zero at the wall to 
some supersonic value a short distance from the wall. A weak wave is 
incident upon the boundary layer and the reflected waves and upstream 
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influence were evaluated. Now with the improved model, however, both 
investigators found, contrary to experimental evidence, that the upstream 
influence is negligible. Although Robinson was chiefly concerned with 
the determination of the upstream influence, Lighthill considered, in 
some detail, the local reflection of the shock wave. He concluded that 
a shock is reflected locally as a "pressure ridge," that is, a rapid 
compression followed immediately by a rapid expansion. Now a note was 
added in the proof of reference 8 to the effect that experimental results 
of Mair and Bardsley (reference 9 ) } obtained at the Fluid Motion Labora- 
tory, Manchester, establish that the conclusions of negligible upstream 
influence and the shock being locally reflected as a pressure ridge are 
correct for the reflection of a weak shock from a turbulent boundary 
layer. (These results, however, had been previously observed by Ackeret 
and Liepmann.) As Lighthill pointed out, the conjecture that the theory 
is correct for velocity profiles typical of turbulent boundary layers 
but incorrect for profiles typical of laminar boundary layers may be 
valid. But until such time as the theory is modified to account for vis- 
cosity, as suggested by Lighthill, the present theory is incapable of 
predicting the effects of the interaction between a shock wave and a 
laminar boundary layer. 

The preceding theories have neglected viscosity and heat conduction 
while considering infinitesimal waves and small disturbances. It would 
appear that the assumption of infinitesimal waves and small disturbances 
is valid since the linear theory, to be a good approximation, requires 
that the slopes of the streamlines be small. For weak incident shock 
waves, the slope of the displacement thickness would probably be small, 
except perhaps in the i m mediate vicinity of the wave. The neglection of 
viscosity, on the other hand, seems to be more serious. If the effects 
of viscosity were not too important, one would expect that, for fixed 
Mach number and shock strength, the effect of Reynolds number should be 
rather small. Actually, this is not the case since the results of Barry, 
Neumann, and Shapiro clearly show that changes in Reynolds number have a 
marked effect upon the shock-wave - laminar-boundary-layer interaction 
(reference . Hence it appears that a theory capable of predicting the 
effects of this complicated phenomenon must include the effects of 
viscosity. 

At present, as far as the authors are aware, the only solution that 
includes viscous effects was given by Lees (reference 10) . The Lees 
theory is based on an approximate relationship between the slope of the 
displacement thickness and the external pressure gradient derived from 
Prandtl-Meyer flow. By combining this relation with the momentum- 
integral equation and neglecting terms of the order of Re - -*- Lees arrived 
at a third-order, linear, differential equation for the pressure. On the 
assumption that this equation is valid all along the boundary layer, he 
was able to determine the boundary-layer growth and pressure distribution. 
He found the upstream influence comparable with that which is observed 
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experimentally and showed, also, that the boundary layer always sepa- 
rates except in those cases where the incident wave is quite weak. How- 
ever, his theoretical pressure distributions fail to exhibit the charac- 
teristic downstream behavior. This point will be discussed later. 

Following the approach of Oswatitch and Wieghardt in their study of 
the growth of disturbances in a supersonic stream outside laminar or 
turbulent boundary layers (reference 11), the problem is treated as an 
outer flow, with a shock, in equilibrium with a boundary-layer flow. 
Instead of expressing the deflection of the streamline in terms of pres- 
sure rise, as in Lees 1 problem, a procedure after Oswatitch and Wieghardt 
is taken by connecting the vertical velocity of the outer flow with the 
streamline deflection, as given by the momentum integral of the boundary- 
layer flow. This reduces the problem to an inviscid one which can be 
solved systematically to any order of approximation. 

This study was conducted at the Graduate School of Aeronautical 
Engineering of Cornell University under the Sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 


SYMBOIS 


a speed of sound 

A,B,C,D,E constants 

b = &!*/ 6 0 * 

Cp, c v specific heats at constant pressure and constant volume, 

respectively 

f,g functions of i - m^H and £ + m^Ti, respectively 

F = 2t| - 2t|3 + rfr 

defined by equations (8) 

g = |nd - n ) 3 

h function of I - m^T} 

H = b*/e 


I 
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k 

K 


coefficient of thermal conductivity 
✓ / 

Karman momentum equation for compressible viscous fluids 


m = ]/ m 2 - 1 

M 

Mo = u/a Q 

M co = u/a M 

P 

P' 

Pr 


R 

Re 


Re x - U x /\»o 
T 

u,v 


u',v* 

U 


x,y 


Mach number (u/a) 


pressure 

perturbed pressure 
Prandtl numb er (pcp/k) 
universal gas constant 
Reynolds number (u5 0 */Vq) 

temperature 

velocity components parallel and normal to flow direction, 
respectively 

velocity perturbations nondimensionalized with U 
free -stream velocity 

coordinates parallel and normal to flow direction, 
respectively 





* 
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8 

8 *- 


8o* 


8' 
8 * * 
A 


ratio of specific heats (cp/cv) 

boundary-layer thickness 

boundary-layer displacement thickness 

undisturbed boundary-layer displacement thickness 

boundary-layer thickness in transformed plane 

boundary -layer displacement thickness in transformed plane 

disturbed boundary-layer displacement thickness nondimen- 
sionalized with S Q * 


e flow-deflection angle 

r\ = y/8 r or y/5 0 * 

boundary-layer momentum thickness 


9 

9' 

Xi, \2, X 3 


[l* i 1 ■ £ 


\_pu_ dyj 
Ue/Peiie / 


boundary- layer momentum thickness in transformed plane 
roots of equation ( 16 ) 




^ U e 

dx 


n 

V 


I = x/8 0 « 


coefficient of viscosity 
kinematic coefficient of viscosity 

density 


a = 1 - 


7 - 


- Mo 2 ' = (l + J ~~- I^ 2 ) 1 


shear stress 

perturbation velocity potential, nondimens ionali zed 
with U8 0 * 


0 
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$ total potential 

X function of x f and y* 

^ stream function 

Subscripts 

d evaluated at distance of upstream influence 

e evaluated along edge of boundary layer 

1 1, 2, or 3 

2 laminar case 

o standard condition, such as stagnation point or undis- 

turbed state 

s separation point 

t turbulent case 

v evaluated at wall 

1,2,3 evaluated in regions 1, 2, or 3^ respectively (see fig- 1) 

oo evaluated in undisturbed free stream 

OUTLINE OF PRESENT INVESTIGATION 


In the present study the effect of an infinitesimal compression 
wave incident upon the laminar boundary layer along a flat plate in 
supersonic flow will be considered. The pressure gradient in the flow 
direction is determined to the order of approximation of the ordinary 
boundary-layer theory solely by the 'shape of the boundary. The poten- 
tial flow thus completely determines the boundary-layer growth. In the 
problem of interaction between a shock wave and a boundary layer, how- 
ever, the phenomena are quite different. In that case the pressure dis- 
turbances that are propagated upstream in the subsonic portion of the 
boundary layer will affect the boundary-layer growth. It becomes clear 
that the outer flow influences the boundary -layer growth and that the 
boundary-layer growth, in turn, influences the outer flow, so that a 
solution that is simultaneously compatible with the two flows must be 
sought . 
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Another significant difference between the interaction and noninter- 
action problems occurs in the separation phenomenon. For the usual case 
of boundary-layer flow against a pressure gradient, it is known that the 

flow separates from the wall at the point where = 0. Downstream 

Ww 

of the separation point the original flow is deflected away from the wall 
by the backflow, a vortex layer separating the two flows. The interac- 
tion of a shock wave with a laminar boundary layer, however, produces a 
different form of separation. Experimental evidence indicates that along 
the wall there exists a short thin region, extend in g slightly upstream 
from the foot of the shock wave, in which the flow has practically zero 
velocity; the boundary between this separated flow and the main boundary- 
layer flow is a vortex sheet. A plausible explanation of this phenomenon 
is proposed by Lighthill, who suggests that the pressure discontinuity 
occurring when a shock interacts with a laminar boundary layer causes the 
flow to separate into some form of "bubble" at the base of the shock wave. 
Once such a bubble is produced, the boundary layer upstream of the sepa- 
rated region will be deflected so as to increase the external pressure 
gradient. This results in further separation of the boundary layer and 
hence further upstream deflection. This process of repeated separation 
causes the edge of the bubble to move upstream until such time as the 
induced pressure gradient is no longer able to cause separation of the 
boundary layer. The net result is the one indicated in experiments. 
Immediately behind the point of interaction, the deflected flow, in 
returning to its original direction parallel to the wall, undergoes a 
very rapid compression, thereby causing an extremely large pressure gra- 
dient. It is known from experimental results that the boundary layer 
downstream is turbulent, but whether the extreme pressure gradient causes 
separation and hence causes the boundary layer to become turbulent or 
whether the flow at that point is already turbulent and thus can sustain 
large pressure gradients must remain at this time a matter of conjecture. 
Whatever the case may be, in the problem for weak shock the flow will be 
assumed to remain laminar downstream of the point of incidence (cf. ref- 
erence 12) . Positive infinity in relation to the problem will be taken 
to mean a distance downstream of the order of 100 boundary- layer thick- 
nesses. In addition, any mention of boundary-layer separation will refer 
to the "dead air bubble" phenomenon. 


In the present theory it is assumed that the boundary-layer equa- 
tions are applicable over the whole plane. Strictly speaking, the valid- 
ity of the usual boundary-layer assumptions near the base of the shock 
wave is doubted, since separation- occurs. However, it is believed that 
for weak incident waveB, a theory based on the above assumptions will 
still enable one to determine some of the important characteristics of 
the flow. In particular it is desired to determine the extent of the 
separated region, since this investigation is merely the first step of a 
proposed theory which will endeavor to take separation into account. 
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The more refined theory is then expected to furnish, in the separation 
zone, slight corrections to the simplified theory. 

The flow field is assumed to he divided into two distinct regions 
as follows: 

(a) A thin layer adjacent to the plate wherein there is laminar- 
type houndary-layer flow 

(h) An outer supersonic potential flow field 

Since an infinitesimal compression wave is assumed, the upstream 
pressure propagation is expected to perturb the boundary layer just 
slightly, thus imposing small disturbances upon the outer flow. This 
justifies the linearization of the supersonic outer flow. (Results of 
this theory will of course be applicable only for fairly weak shock 
waves.) The flow field has been assumed to be separated into two dis- 
tinct regions. In reality a continuous variation in velocity from zero 
at the boundary to some main-stream value a short distance from the wall 
is known to exist. The actual distance at which the velocity reaches its 
main-stream value is rather indefinite; hence the so-called boundary- 
layer displacement thickness S*(x) is defined. Physically it represents 
how much the potential flow streamlines are deflected because of the 
reduction in mass flow caused by the retarded velocities near the wall. 

The viscous effects are thus characterized by a layer whose thickness 
is 8*(x), and the link connecting the outer flow to the boundary-layer 
flow is the condition that the line y = 8*(x) is a streamline of the 
outer flow. The treatment of the boundary layer in the present study is 
based upon a momentum-integral eguation previously derived by Howarth 
(reference 13). The deflection of the boundary layer is expressed in 
terms of the vertical, velocity component along the edge of the layer and 
this relation is used as a boundary condition for the outer flow. Hence 
the over-all problem is then reduced to the solution of a purely super- 
sonic potential-flow problem. 


LAMINAR -BOUNDARY -LAYER EQUATIONS 


The Prandtl boundary-layer equations for the two-dimensional steady 
flow of a viscous compressible fluid are: 


ox dy p 8x P 8y\^ 8y/ 
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0 


1 8p 

P by 


~— ( pu.) + ^-(pv) = 0 
ox ay 



p = RpT 


Here u and v denote, respectively, the x and y components of the 
velocity and p, p, and T, the pressure, density, and temperature of 
the fluid. In order to simplify the problem, Howarth (reference 13) 
assumes that the coefficient of viscosity p. is proportional to the 
temperature and the Prandtl number is unity. These assumptions, of 
course, introduce some approximation, but according to the results of 
Emmons and Brainerd (references lh and 15), the boundary- layer character- 
istics depend very insensitively on the form of p(T) and the change of 
Prandtl number for moderate Mach numbers. Therefore, the simplifications 
brought about by this assumption justify the slight Iobs in accuracy. 

By introduction of the stream function i|r 


Pu = P 0 ^- 

oy 


PT = _ Po 21 
ox 


and by introduction of the transformation 


y* = 


X 1 = X 





* = 



X(x’,y') 
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the momentum equation is transformed into 


d 2 X dX 
dx* by 1 by* 


a 2 x ax 

(by*) 2 ax« 




_JL_ X - d 2 - X ■ + v n -~ 3 - 

2ae 2 (dy 1 ) 2 (8y')3 


where p Q , T Q , and v Q stand for the density, temperature, and kinematic 

viscosity of a standard condition, the stagnation point, say, and are 
constants, and the subscript e indicates the condition at the edge of 
the boundary layer. In the case of a thermally insulated plate 

= 0, the condition of uni ty Prandtl number admits a special solu- 
tion on the plate (reference 16 ) 



c pT + | u 2 


CpT e 



Upon elimination of T, there results the equation: 

d 2 x ax a 2 x ax a 3 x i 

8x' 8y' by' (dy 1 ) 2 8x' ° (8y*) 3 


^e 


du e 

dx’ 


r • 

_ 

« 

k . 7 - 1 

2 fbx \ 2 

^ " W 

7 8 2 x ^ 

[ «s 8 

2a e 2 (dy') 2 
(1 


dUo 

Thus, it is seen that when — — = 0, that is, at constant pressure, this 
equation is identical to the equation for the incompressible fluid. 

Let 6' be the boundary-layer thickness in the transformed plane, 
and define the corresponding displacement and momentum thicknesses by 
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Then, by integrating equation (l) across the boundary layer, tbe equa- 
tion relating tbe boundary-layer thicknesses with the velocity is 


ae* 

dx' 


+ 






1_ 

u e dx' 



( 2 ) 


The similarity of this equation to 'that of the incompressible flow sug- 
gests the applicability of the Karman-Pohlhausen procedure. 

If the velocity profile is given by (reference 13) 

= F(n) + AG(tj) 


where 


with 


it follows that 


F = 2tj 


2t]3 + 


G = | n (l - n ) 3 





0* 




S'* = 36 

- A) 

1 


V 

(3) 

II 

: + t) 
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Substitution of the quantities of equations (3) into equation (2) 
gives rise to a relation between the boundary-layer thickness 6' and 
the velocity Ue( x) , with a parameter A. Now, in the present problem, 
as the flow in the boundary layer must be compatible with the potential 
field, the flow deflection dS*/dx of the boundary layer determines the 
flow in the outer field. Therefore, 8' should be eliminated. For this 
purpose, let the displacement thickness be introduced: 



But in the transformed plane there is 



By means of the temperature -velocity relation and Pohl hausen 1 s profile 
this will yield 


5‘* 


© 


1/2 


T 

_£ 

T e 


5* + 


Me 


U + 361 + 1 _1L_ 

V 630 3 2520 


A + 



Now, from the second of equations (3)> S'* is already expressed in 

terms of 8* and A. The elimination of S'* , then gives 


7 - 1 Me 2 
72 252 


_L du e / 

V 0 dx \ X 


7- - 1 
2 



(S ') 5 


+ 


JL dUe fi 7 - 1 M 2 \ ( 7 - 1 71Me 2 1 V n ; 

V D dx V 1 2 ^ JV 6 2 520 12oJ^ S ' 


7 - 1 367Me 2 
. 2 630 



w 
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Theoretically, this equation defines 8' as a function of 5*. 
Only in the case of small pressure gradient, namely, 

Ue = U(l + u* + . . . ) 


u* « 1 


however, is a simple analytical solution possible. For by retaining 
only linear terms in u* and its derivatives there is obtained (refer- 
ence 17) 


U /l _71_ ff -2 
v o \3 2520 




+ 


37 

315 


263 

630a 


263 

630 


(7 - l)^ 2 

a2 



6' = 


2-7 

a 2(7-l) _ 2-7 
2 


4-37 

ij 


U‘ 


8* 


( 5 ) 


where Mq = 


U 


and 0 = 1 


Mo . Moreover, since u' is assumed 


to be small, that is, for weak interaction, the boundary- layer thick- 
ness 8' must also differ from the unperturbed value by a small amount. 
Therefore, to the first-order approximation there results the solution: 


S' 


2-7 

02 ( 7 - 1 ) 


- _3I 
315 


+ 


2 ^3 \ _ 

630a) 


263 (7 - i)Mq 2 /263 
630 6-57 U3O0 





4-37 

+ 2^_7 « 2^(7-l)/263_ 
2 ^6300 



+ 


\ 1 3( 2-7) 

JL/l Tier" 2 _ a^\/263_ _ jyC\ h ^T) Ml ( 5 *)3 
v o\3 2520 945/V63O0 315/ dx 


( 6 ) 


Let 8 0 * be the unperturbed displacement thickness and write 


8* = S Q *(1 + A) 
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Substituting equation (6) in the momentum eguation (2), replacing x* 
by x, and keeping only first-order terms, one obtains, by taking into 
dS* 

account 5 0 * __ = Constant, the following lihear equation (refer- 
ence 17): 


where 


1 _ d_A 
gp d| 



+ 


3g 3 


du* 

d| 


+ 5gij. 



(7) 





82 - -81 




g 3 


i _ 7b ff 2(7 - l) 
263 

(o _ 7)M 2 l^ + JZ^88l 

= U 7jM o _ 2 9731 3600^7 P 

263 


2a 


(8) 


g , . -b e a7=i-L.f Sorfi - JCt a 

4 e 888l263j l 263 


-3 


with 
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The momentum integral thus relates the displacement-thickness perturba- 
tion A to the velocity perturbation and its derivatives along the edge 
of the boundary layer. 


SOLUTION OF BOUNDARY- VALUE IROBIEM OF OUTER FLOW 


Until now the study has been directed at obtaining an explicit rela- 
tionship for the boundary-layer displacement thickness as given in equa- 
tion (7)- Now the extent of the interaction zone is only a fraction of 
the distance from the leading edge of the plate, hence 5 Q *(x) can be 

assumed to be constant over this range. Equation (7) then becomes a 
linear equation with constant coefficients. 

The problem is now one of determining the supersonic potential flow 
past a thin body whose thickness at any point corresponds to the perturbed 
boundary- layer thickness at that point. Consider the problem in the 
coordinates shown in figure 2. 

The undisturbed thickness is seen to correspond to y = 0. Strictly 
speaking, the undisturbed thickness should be at y = 5 0 * in the physical 
plane. But since the outer flow does not depend upon the location of the 
origin of coordinates, the boundary conditions can be simplified slightly 
by choosing the axes as indicated in figure 2. Then, consistent with 
linear theory, the boundary conditions are satisfied at the undisturbed 
surface, namely, y = 0. 

The linearized differential equation for the perturbation velocity 
potential of the supersonic stream is 


nW 


2 d g 0 
St 2 


- = o 

8t] 2 


(9) 


where 


i = yr^rr and 
6 o* 


7 = 


_ JL _ 


V* 


and the general solution Is given by 


0 = fU - =Wl) + g( I + nwi) 


( 10 ) 
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The incident waves will be taken as 


Sit + IDwTl) = 0 

gU + mjl) = --|-(i + m^) 


when 

when 


£ + m^T] < 0 
| + m^T) > 0 


( 11 ) 


This corresponds to a simple compression wave, of deflection angle e, 
incident upon the origin. (According to the above definition, the abso- 
lute value of e will be used in any computations.) 


Upstream Solution 

The incident wave causes the boundary layer to be perturbed upstream 
and downstream of the point of incidence. Since no disturbances can be 
propagated upstream into the supersonic flow, the physically possible 
solutions in region 1 will be waves of the form = f ( £ - m^T]) . The 
normal velocity component at the boundary must be zero, hence the condi- 
tion that the slope of the streamline is equal to the slope of the dis- 
placement boundary-layer thickness 5* is imposed. The velocity vector 

of the supersonic flow has 1 components [u(l + u' ),UvQ, where u' = 
and v' = c>0/dT}, so that the linearized boundary condition becomes 


or 


d A = 

di 


at ^ = 0 


= -nicof’U) at ri = 0 


( 12 ) 


Now since the boundary condition is on d A/d| rather than A, 
equation (7) is differentiated with respect to |, considering 6 0 * to 
be constant: 
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a a . a 2 A 


du’ 


d 2 u' 


d3u’ 


gl Sr + iis - Ir + g 3 sir + jr 


1 d 2 A t , . „ du' . = d^u' 

- ^(g 2 u. + 3g 3 jp ♦ 5*4 ^ 


1 d A ( du' , _ d^u' , _ d^u 

Si (f 2 a£ 3g 3 ai 2 564 d^3 


d^uA 
a?3 ; 


(13) 


If the "boundary condition (12) is substituted into equation (13), and it 

is noted that at tj = 0, = f'*(£), and so forth, an ordinary 

d£ 

nonlinear differential equation for f(|) will be obtained. By line- 
arizing this equation: 


Sk f±V (t) + g 3 f , “U) + (g 2 + m a,)f * *( l) + gim^f^l) = 0 (14) 


The coefficients in equation (l4) are constant for a specific case of 
Mach number and Reynolds number, hence the solution can be immediately 
written as 


fU) = Ao + Ae 7 ^ 1 + B 0 e A l^ + C D e A 3 & 

where the Aq's are roots of the equation 

S^ 3 + g 3 A 2 + (g 2 + m«,)A + gjm^ = 0 


(15) 


(16) 


For the cases under consideration equation (l6) has three real roots, 
one positive and two negative. Let A 2 > 0, Aj < 0, and A 3 < 0. In 

order to have' all disturbances vanish as | — > -» it is seen from 
equation (15) that 


B o = C o = 


= 0 
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The upstream solution is then 


f(S- - mji) = Ae 2V 


(17) 


or in the physical plane 


* - J^ y) 

x - m^y < 0 


(18) 


The solution in region 2 follows immediately. The equation for the outer 
flow is linear, hence, by adding solutions, 


02 = Ae ° 


e 1 


®eo 5 o* 


x + m^y) 


x - m^y < 0 
x + m^y > 0 




(19) 


Downstream Solution 

In region 3 the solution to the linearized supersonic flow equa- 
tion is 


03 = g( i + nwi) + h( I - m ro Ti) 


But g( | + MooT)) is known since it is the incident wave, hence use of the 
boundary condition at the interface will result in an ordinary differen- 
tial equation for h(t). The nondimens ionalized boundary condition takes 
the form 


fr = - m o£'(^) - € 


at t) = 0 


( 20 ) 
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since 


Si] 


= -e - 


Substituting equation (20) into equation (13)., 


gl 4 h lv + g 3 h* r ' + (g 2 + noo)^' + gpnUa* = Tg x e 


( 21 ) 


It can be seen that the complementary solution will be identical to 
equation (15) • Hence the complete solution is 


Ap £ Aq_ ^ Ao I 

h(£) = D 0 + Ee + Be + Ce J - — I 


( 22 ) 


where again A2 >0, Ai < 0, and A3 < 0. Choose D 0 = 0 without 
affecting any of the physical quantities such as velocity or pressure 
since h(|) merely represents a velocity potential. To eliminate the 

possibility of any velocity becoming infinite as £ ► a> put E = 0. 

Then the solution becomes 


A-j£ Ao£ p 

hU) = Be X + Ce J ~ I 


, ^^(^-"Wl) „ ^3(6-^) g 

h(l - m^Ti) = Be + Ce - 5 - 


mooli) 


(23) 


Then the perturbed velocity potential in region 3 is 


0 3 = Be 8 ° + Ce 5 ° 


2e _x_ 
n\» 5 o* 


x - m^y > 0 


L (24) 


Determination of Constants 

The problem now involves four constants A , B, C, and D^ with 
three parameters M^y g, and Re. Additional conditions must therefore 
be imposed in order to determine the constants as functions of the 


•kche constant D is a measure of the downstream boundary-layer 
thickness as shown in the following discussion. 
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parameters. It is noted that the outer flow field was divided into dis- 
tinct regions by the lines OS and OM (figure 2), so that the boundary 
layers in regions 1 and 3 were treated independently of each other. 

These boundary layers are therefore related at the point of incidence by 
the following conditions: 

(a) The boundary -layer displacement thickness must be continuous 

at x = 0 

(b) The wall pressure must be continuous at x = 0 

px o 

(c) Discontinuities, if any, must satisfy lim / K dx = 0 

Xq * 0 J -Xq 

where K = 0 denotes the Karman momentum equation for compres- 
sible viscous fluids 

Condition (a) .- Compute the boundary layers for regions 1 and 3 and 
then match the displacement thicknesses at the origin. One has as a 
boundary condition 


at t) = 0, i = 1 or 3 


is the total boundary- layer displacement 


From equation (l8) one obtains for region 1 

d Ag £ 

= -m^Ae 

A2I 

= -m^Ae + Constant 

All disturbances vanish far upstream so that Ap(-oo) = 0. But since 
>2 > 0, the exponential term vanishes as i — > - 00 , hence the constant 
equals zero . Therefore 


d A^ 80^ 
d£ 


8j* 

where A± = g— ^ - 1 and Si* 
thickness in regions 1 or 3- 
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For region 3 one obtains, according to equations {2k ) , 

= -m^BAxe 7 ' 1 ^ - mj^e^ 1 

At | Aq | 

A3 = -m^e x - m^pe J + Constant 

Let ^3(00) = D and since the exponentials vanish at positive infinity, 

the constant equals D. 

Therefore 


or 


A3 = -m^e 7 ^ - m ro Ce A 3^ + p 


Ai A3 

83* = 8 0 *[1 - m^e 0 - mjle 0 + D| 


Thus it is seen that for the total downstream thickness 


i 



= B 0 *(l + D) 


By putting Sp* = 83* at 


x = 0, 


A = B + C 



(25) 


Condition (b) . - Since the validity of the usual boundary-layer 
assumptions in the interaction zone has been assumed in the case of weak 
shock interaction, the pressure at the wall, as a first approximation, 
equals the pressure at the interface. The pressure at the wall as well 
as at the interface is therefore required to be continuous. As a conse- 
quehce, the pressure jump due to the incident wave must necessarily be 
neutralized by a reflected expansion wave. This deduction is act uall y 
confirmed by experiments. 

The pressure in each region consists of an undisturbed pressure p 0 
plus a perturbed pressure p' due to the incident wave; P = P Q + P l • 

Since the undisturbed pressure is the same throughout the flow field, 
only the perturbed pressures, which, by small perturbation theory. 
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are p' = -PooQu’, will be considered- The velocity potentials 

and have previously been found, and, accordingly, the pressures are 

given by 


Pi' = -iM^QA^e 


P3 f = -rM^o 


^1 


BAt e 


u . n -v _ U 


+ CA ^ 6 


2e 

®oo 


The condition for continuous pressure at the origin is 


(26) 


(27) 


Pi' = P 3 ‘ 

hence 


AA 2 = BAi + CA3 - |i (28) 

This condition shows that the discontinuity between regions 2 and 3 cor- 
responds, in linear theory, to an expansion of the flow from regions 2 
to 3 By a wave of the same magnitude as the incident wave. Now since Ap 

and A3 &re negative, the downstream pressure along the interface y = 0 
is 


p 3' 




m„ 


P 0 (2e) 


(29) 


as £ — > «>. This asymptotic value of the pressure is exactly twice the 
incident pressure rise; therefore, as anticipated, the incident wave is 
reflected as a regular reflection from a solid boundary. 

Condition (c) . - The Karman momentum- integral equation for compres- 
sible fluids is 



pu c 


ay - ue H 



dp 

pu dy =-^ S 


'v 
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where t w indicates the wall shear stress. With the ordinary definitions 
of the displacement and momentum thicknesses, the equation becomes 

^(Pe^e 20 ) + PeUeS* ^ = t v 


From the outer flow 


du, 




e e dx 


dp 

dx 


thus 


dx 


(p e u e 2 e) - 8* 


- T„ = 0 
dx w 


From the previous notation it is seen that 


K = A(p e u e 2 0 - 5*p) + p 


d5* _ 
dx~ w 


As the boundary-layer regions 1 and 3 are independently considered, 
the pressure gradient and velocity therein might be discontinuous at the 
origin. In order that the discontinuities be consistent with the dynami- 
cal equations, it is necessary that: 


or 


lim 

Xq ^ 





lim 

Xq >>0 


^PeUe 2 © - 5*p^ ° + f p dx“ dX "/' Tvdx 


= 0 


Now 5* is required to be continuous at the origin, 
but finite, and the shear stress t w , finite; then. 


■^t, discontinuous 
dk 

Kh t=» H m-f 
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t wo integrals vanish and 

lim ( PgUe^S - 8*p\ ° =0 

x o >0\ J-x Q 

remains. Therefore, in the limit, 

P^u^ 2 ^ - PjU-^fi^ - S*(o)^P2 - = 0 


(30) 


The condition that = p^ at x = y = 0 has already been imposed 
and hence equation (30) reduces to 


Pl u l^®l = P 3 U 3^®3 at x = y = 0 


(31) 


This equation, linearized in u* and its derivatives, becomes 


a A*(lr ' 1 &) = °e( u 3' - V) 


where 


a l = “ 


Re 


/ 37 263 \ 3 

\315 " ^30?/ 


3(2-7) 

,2(7-1) 


Ills 


i/JiV 2 . sf\ 

_3V2520 r 955] 


( 37 263 \ 

\315 6300; 


(32) 


2-7 


a, 3 - f£_\ 

2 \315 630 0/ 


2 + - x ) 


The pressures were equal at the origin, so that, in linear theory, 
streamwise perturbation velocities are also equal. Therefore Uj_' 

at x = y = 0. 


the 
= u 3 « 
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Now since £ 0 and cu> ( 0, for tie range of and R e 

studied 


dU}_ 1 du^ ' 
dx - dx 


at x = y = 0 


(33) 


It is therefore seen that, in the linear theory, continuous pressure 
also implies cont inu ous pressure gradient. Since 

duj* A/v^ 2 

"to" = 57* at x = y = 0 

duo ' 1/2 2\ 

UbT - SJv’' 1 + ) at x - y • 0 

there is obtained from equation (33) 

AA 2 2 = BA X 2 + CA3 2 (3^) 

Conditions (a), (b), and (c) lead to equations (25), (28), and (3^-) 
for four constants A, B, C, and D. But previously it has been shown 
that 

lim 83* = 8 q *(1 + D) 

| >• 00 


Thus it is seen that D merely determines the downstream boundary- layer 
thickness, whereas A, B, and C determine the local character of the 
perturbation. Therefore, if the downstream thickness can be estimated, 
there will be a deter mina te problem. Assume for the present at least, 
that D is determined. The constants A, B, and C can then be solved 
for" from the following equations: 

to 2 = B *1 + - If 


A = B + C - 


m_ 


AA2 2 = BT^ 2 + CA3 2 
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The. solution is 


A 


B 


C 


~( A 3 + A l) / 2e\ + A 1 A 3 (_ D_\ 

( A 2 “ ^ 3 ) (^2 " A l)l ( A 2 “ A 3)( A 2 " A l) V ™°°/ 

(A3 + A 2 ) / 2e\ V ^3 /_ D_\ 

(a 3 - a 1 )(a 2 - Aj \ “J " (a 3 - A-l)(A 2 - A x ) V “W 

-(a 2 + Aj) / 2g \ + A 1 A 2 / D \ 

(A3 " Ai)(A 2 - A 3 ) V m oo/ (A 3 - Ai)(A 2 “ A 3 ) V m coj 


-'v 




(35) 


Evaluation of D .- It is seen from equations (35) that the local 
perturbations are determined once the constant D is specified. As 
it has been shown that D characterizes the thickening of the boundary 
layer through a shock, it can be determined approximately by the 
following consideration. 


When a weak shock is incident upon the boundary layer, it is assumed 
that the Karman momentum- integral equation 


d0 

—■ + S 
dx 


_1_ 

+ 1 

2 -4- §2 

I 1 

dUg 

Pe 

dx \ 

0/ 

u e 

dx 


w 


P e u e 


is valid throughout the disturbed boundary-layer flow. For a given pres- 
sure distribution, the growth of the boundary layer is governed by this 
equation. When a shock is incident upon the boundary layer, it has been 
theoretically predicted and experimentally verified that the flow up- 
stream of the point of incidence is separated over a considerable portion 
of the disturbed flow; hence the shear stress becomes relatively unimpor- 
tant. Downstream of the point of incidence, since the flow is dominated 
by pressure forces, whether the flow is laminar or turbulent the shear 
is known, from experimental results, to be very small (reference 18 ) . 

Hence for the present approximation, the shear stress can be neglected. 

Furthermore, if the ratio 6*/0 = H is regarded as a parameter, 
the momentum equation gives the momentum thickness as a function of u e (x) . 

Now, it is known that H is substantially increased by the shock, and 
to a lesser extent by the compressibility. In going through a shock, H 
first increases because of the shock and then decreases because of the 
drop in Mach number. Hence, as a first approximation, H can be taken 
to be some constant average value over the interaction range. (This is 
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similar to the procedure mentioned, by Nitzberg and Crandall (refer- 
ence 19)-) In particular, H will be taken as the value ahead of the 
interaction. 

The momentum- integral equation now becomes: 


1 d0 + _1_ 

0 dx p e dx 


(2 + H)— 

u^ dx 


= 0 


or 

~^log e 0p e Ue 2+H ) = 0 

Integrating between regions 1 and stations far upstream and downstream 
of the point of incidence. 


2+H ^ 2+H 

e l p l u l " e 3 p 3 u 3 


or 


!£ fi 

8 1* p 3\ u 3/ 


(36) 


For laminar flow over a flat plate where p « T, P r = 1, 

and (8T/8y) w = 0, the formula given by Lees (reference 20, p. 119) 

for H can be taken: 


H = 2.50 + 3-50 



The thickening predicted by equation ( 36 ) is shown in figure 3* 
Values of the density and velocity ratios were computed exactly and to 
first order in the deflection angle. For small deflection angles, the 
agreement between the "exact" thickness ratio and that obtained from 
linear theory is quite good. The crosses on figure 3 are the experi- 
mental values for various Reynolds numbers obtained by Barry, Shapiro, 
and Neumann (reference k), and there is fair agreement between theory 
and experiment. Actually, the experimental data are visual est ima tes of 
the ratio of boundary-layer thicknesses, so that for a comparison between 
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theory and experiment, the ratio of houndary-layer thickness to displace- 
ment thickness would have to he assumed to he the same upstream as well 
as downstream of the point of incidence. In view of the assumption and 
the allowable experimental error in the visualization of schlieren photo- 
graphs, the experimental values must he regarded as qualitative. The 
theoretical values of the thickness ratio are near the lower limit of the 
experimental values. This is prohahly due to the fact that the value 
of H was underestimated, since an average H through the interaction 
range would he larger than the initial value. 

Evaluating D according to the linear theory, 


and 


hence 


But 


Therefore 


s 3* = fl/M 
s l* p 3\ u 3/ 


Pi 

p 3 


1 




— = 1 + — 

U 3 




1 + D 


D 


2e4 

m \ 
00 ' 


2 + H 
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Now since 


H = 2.50 + 3-5o( 2 ^-^)m co 2 

D = “ °' 3M “ 2 ) (37) 

m oo 

Effects of Downstream Thickening 
Consider the constant A governing the upstream flow: 

A = ~(*3 + Ax) / 2 e\ + *1*3 ( D \ 

(Ag - A 3 )(A 2 - 7 , 2 ) ' ~ ^ 3)^2 “ *i)\ m °°' 


or 


A = - — 


2e 

®00 


- 

1*3 

+ *1, 

*1*3 

(k.50 - 0.31^ 

(a 2 - 

*3) 

(*2 ‘ 

*1) (*2 ~ *3' 

(*2 " A l] 

V / 


The coefficients are given in table 1} hence it is seen that for all 
cases in the present range of Mach nuiribers and Reynolds numbers the con- 
tribution of D toward A amounts to less than 10 percent. Thus it is 
seen that, in the present range, a precise knowledge of the downstream 
thickness is not essential to the determination of the upstream flow . 

The effect on the downstream flow is not quite so straightforward, 
since two constants B and C are involved. As an indication of the 
downstream flow, consider the pressure distribution along the wall. It 
has been seen that, far downstream, the pressure attains a constant 
value, namely the pressure that would have been anticipated had there 
been a regular reflection. Immediately behind the point of incidence the 
pressure is much lower than this end pressure since the flow has just 
undergone an expansion to the pressure that existed before the point of 
incidence. The manner in which the pressure proceeds from its value at 
the origin to its final downstream value is then of some importance. 

Does the pressure increase monotonically to its final value, or does the 
pressure at any point rise higher than the final value? Liepmann (refer- 
ence 3) has found experimentally that for the reflection of a shock wave 
from a laminar boundary layer there indeed exists a definite downstream 
overcompression. The pressure rises sharply to some value higher than 
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the end pressure and then tapers off. The experimental findings of Barry, 
Shapiro, and Neumann (reference 4) also confirm this observation, although 
their results show the overcompression to be slightly less pronounced. 


• Consider the behavior of the theoretical downstream pressure distri- 
bution. Along the interface the total pressure is given as 


p 3 = p Q - (sApe 1§ ° _ + CX 3 e 3 S o _ |e 

“ooy 


(38) 


and 


, 2 7 M^ 

p i + 


(39) 


as 



I — 


so that far downstream, consistent with linear theory, the correct pres- 
sure behind the regular reflection is obtained. Then 



wherein e > 0 and D > 0. Now A2 >0, Ap < 0, and A 3 < 0, where A2 
and Ap are of the same order of magnitude, while A 3 is about ten 
times as small. Hence B < 0. 

Consider 

c = -(A g + Ap) / 2g \ + A 1^2 / D \ 

(A3 - Ap)(A 2 " A3) \ “«/ (A 3 - Ap) ^A 2 " A3) \ 


26 

1 

[A 2 + 

A x ) A l A 2 

/4.50 - 0.3^^ 

“co 

( A 3 ' 

-M 

(A 2 - A 3 ) (A 3 - Ap) 

(a 2 - A 3 ) 

\ V ) 
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Now the second term is always negative in the Mach number range considered 
whereas the first term may be either positive or negative depending upon 
the sign of A2 + Ap. It turns out that for M ro > 2.4 (approx.) the 

first term is positive, hence C >0. For Mg, < 2.4 (approx.) the first 
term is negative, but 


(^2 + A l) 


ApAg A-.50 - 0.3M J\ 

(A 3 - Ap)(A 2 - A 3 ) 

N. 

(a 3 - Ai)(a 2 - a 3 ) \ 11100 / 


hence C is positive. Thus for the whole range of Mach numbers and 
Reynolds numbers, C >0. 

The qualitative behavior of the downstream pressure distribution 
will now be established. The downstream pressure along the interface is 


P 3 = Po + ^ + 7M ™ 2 ?o ^-BApe^lSo* - CA 3 e A 3 5 o*j 

The first two terms represent the constant end pressure while the third 
actually describes the variation of the pressure. Consider the curve 

y = 

It is known that B < 0, C > 0, Ap < 0, and A 3 < 0 m , hence -BAp < 0 
and -CA^ > 0. The constants B and C are of the same order of magni- 
tude, but Ai is about ten times as large as A3J hence 

|bAi| > |ca 3 | 

Therefore, when x = 0, y is negative and, when x becomes large, y is 
positive but small. The point at which the maximum occurs is given by 

= 1 T ( BA 1 2 

A3 - Ap ° Se ( CA 3 2 

TTl ftV \ J 
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If 


y = 


A 

-BApe 


x 

l&o* 


A 

— CA^e 


x 

3&o* 


is now interpreted as the pressure variation, it is immediately seen that 
the theoretical pressure distribution exhibits the overcompression that 
is observed experimentally . 


RESULTS AND DISCUSSION 


It has been seen that if the downstream thickening is estimated, 
the constants A, B, and C are completely deter min ed as functions of 
the flow parameters. With the outer flow thus determined, the boundary- 
layer growth and pressure distributions can now be computed, and hence 
the effects of Mach number, Reynolds number, and shock strength upon the 
upstream influence and the location of the separation point can be 
studied. (Of course, when the results are compared with experiment, the 
Reynolds numbers must be low enough so that the boundary layer will 
remain laminar, prior to the interaction, in the corresponding experi- 
mental case . ) 


Upstream Influence 

In order to estimate the upstream influence of the interaction, 
define a length in which the pressure on the boundary-layer displacement 
thickness decays to a specified fraction of its amplitude at the origin. 
This length would then be determined solely by the exponent A 2 * For 

example, consider the upstream boundary-layer disturbance: 


A 

5 X * = -n \*A8 0 *e 


x 

28 o* 


At the origin 



36 


MCA TR 2868 


If x<i/5 0 * is defined as the distance required in order that 
81* = a^-in^ASo*) at x<j/s 0 * where a < 1, 

A — 

a(-BtoA6 0 *) = -m^AS^e s5 °* 


Therefore 




Since a is a constant it is seen that the upstream influence is then 
inversely proportional to Ap • Values of A2 have been plotted in fig- 
ure 4. (Values of A-^ and A3 are plotted in figures 5 and 6 , respec- 
tively.) For fixed Mach number, x^ S 0 * increases with increasing 

Reynolds number. For large Reynolds number, this dependence can be 
deduced from the equation: 




'82 

— + — A + 

gl). g U gJj. 


The coefficients are functions of Mach number and for the range of Mach 
numbers considered, the roots are real. When the Reynolds number is 
large, it can be shown that 


y^T 

]f~ Zk 


where g^(Re, M^,) < 0 and g^ ocRe. This is in agreement with Lees' 

result. Consequently, for large Reynolds numbers, x^jb 0 * « Re^, For 

fixed Reynolds numbers, the upstream influence decreases with increasing 
Mach number as indicated in figure 4(b). If the disturbance is considered 
to decay to, say, 5 percent of its amplitude at the origin, it is found 
that, at the high Reynolds numbers, and Mach numbers about 2 or less. 


the upstream influence is of the order of 30 boundary-layer displacement 


thicknesses . 
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The result that the upstream influence increases with Reynolds num- 
ber may appear to be disconcerting since the viscous effects are expected 
to be more prevalent at the lower Reynolds numbers and hence to produce 
a greater influence upon the upstream flow. This apparent paradox arises 
because the upstream influence has been measured in multiples of a length 
which is also dependent upon Reynolds number. If the absolute values of 
upstream influence are considered. 


x d 

So* 


<x Re 


1/2 


1 


But 


S 0 * 


Re. 


i 72 


hence 


x d 


« 6 0 * Re 


1/2 


or 


x d 


a Re 


- 1/2 


Thus, measured on an absolute scale, the upstream influence increases 
with decreasing Reynolds number . 

As defined above, the distance of the upstream influence is depend- 
ent solely upon Mach number and Reynolds number. It can be redefined so 
that it will alBO depend upon the shock strength. This is accomplished 
by defining the upstream influence to be that distance at which the dis- 
turbance decays to a given fraction of the undistrubed value. The per- 
turbed boundary-layer displacement thickness ahead of the point of inci- 
dence is given by: 


6 


,* = 


A 


x 

25 0 * 


where 8 0 * is the undisturbed displacement thickness. The upstream 
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influence is then the distance X(j/& 0 * to the point where 5-j* = hS Q * 
and h < 1: 



Now A is negative and decreases with increase of the deflection angle, 
so that for fixed Mach number and Reynolds number, the upstream influ- 
ence decreases with shock strength. If the upstream influence is mea- 
sured to the point at .which the disturbance decays to 5 percent of its 
undisturbed value, b = 0.05. In figure 7 , x^/5 0 is plotted against 
deflection angle. For e < 2°,. the results seem to agree fairly well 
with the experimental values presented in figure 10, reference 2. For 
larger deflection angles, however, the theoretical values are too low in 
comparison with the experimental values. This might be due to the fact 
that the linear theory becomes less accurate as the deflection angle 
increases. 


Boundary -Layer Separation 

Since the outer potential flow is known, the point of separation 
can be computed on the basis that separation occurs when (chi/8y) v = 0. 

Now du/dy is proportional, to du/y ' so that one can just as well use 
(8u/8y') w = 0 as a criterion for separation. In view of the assumption 

of a Pohlhausen velocity profile, in the transformed plane 



hence, as in the incompressible case, separation occurs at A = -12. 

Now it is known from experience in the incompressible case that, in 
regions of retarded flow, the Karman- Pohlhausen method gives values of 
the skin friction that are too high and consequently predicts separation 
too late or not at all. This feature of the method is also to be 
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expected in the compressible case, since it was shown by Howarth that 
the effect of compressibility is equivalent to exaggerating the pressure 
gradient in the incompressible case (reference 13) . In particular, 
Stewartson has used the Karman-Pohihaus en method to treat the case of 
flow against a linear pressure gradient and has shown that the predicted 
distance from the leading edge to the separation point is an overestimate 
of the actual value (reference 21). 


In the present problem, where the flow upstream of the point of 
incidence is subject to a positive pressure gradient, the same overesti- 
mation of the separation point will, of course, be expected. But the 
effects of Mach number, Reynolds number, and shock strength upon the 
separation point are of primary interest, so that inaccuracy in the abso- 
lute values is unimportant. 

It is known from the boundary-layer theory that the parameter that 
defines the separation point is 


A = 


.ism 





By neglecting products of u’ 


and its derivatives. 


A 



(U/crHdu’/dx) 


Wow, from eguation (6), 


to the order of approximation. 


3-27 

Re 5 0 * 7-1 du' 

/ 37 263 \ 2 a dx 

\315 ’ 630 a ) 


Along the edge of the boundary layer, from the upstream solution (see 
equation (l8)) 




du' 

dx 



ho 
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so that 


A = 


^ o 3-27 ^ x 

A Re A 2 7 _ X ^p8 n * 

— a e ^ u 

' 37 263 \ 2 

315 " 630a/ 


For separation 


. ? ' 3i£ 7 . ^s_ 

A Re 7-1 * 2 8 0 * 

- C y / e ^ 0 

/ 37 263 \ 2 

\315 630 a) 


Xg_ 

6 0 * 



/ A Re A 2 ^ 

V 


m 


where 

3-27 

7-1 

R = J 

P — / \ O 

/ 37 _ 263 r 
\315 63O0J 


and A 2 is positive, hence separation ahead of the point of incidence 
will he predicted only when A Re Ap^P > 1* 

Contrary to the separation phenomenon in ordinary houndary-layer 
flow, the location of the separation point is influenced hy the Reynolds 
number. This is easily understood when one considers the fact that 
separation of the flow 1 b controlled hy the pressure gradient of the 
outer flow} according to houndary-layer theory, the outer flow depends 
solely upon the geometry of the hody. Consequently, the location of the 
separation point is independent of Reynolds number. In the present 
problem, however, an outer flow that is compatible with the houndary- 
layer flow has to he found. This relationship is expressed hy the condi- 
tion which requires that the direction of the potential flow he the same 
as the slope of the displacement thickness. Since the houndary-layer 
growth depends upon the Reynolds number, the outer flow, and eventually 
the separation point, must also vary with Reynolds number. 
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In order to determine the effect of shock strength upon the position 
of the separation point, the constant A must he examined: 


A = - — 

“co 


- 

fa 

+ Ai) . A 3 A 3 

A-.50 - 0.3M„ 2 \ 

(A 2 " 

A 3 )^A 2 - A].) (A 2 - A 3 )(A 2 - Ai) 

V ^ J 


For the range of Mach numbers and Reynolds numbers considered, the terms 
in the bracket are positive, since Ag >0, Ap < 0, and A^ < 0. There- 
fore as e is increased A becomes more negative) consequently 
-^AEe A 2 2 f3 becomes more positive and hence the separation point moves 

upstream. Thus an increase of shock strength increases the distance 
between the separation point and the point of incidence . 

To express x s /5 0 * explicitly in terms of Mach number and Reynolds 
number would be difficult; hence this relationship will be presented 
numerically by varying, separately, the Mach number and Reynolds number. 

In figure 8, lines of constant Reynolds number are plotted in the 
xs/So*,€' plane, and the separation point, measured in multiples 
of S 0 *, moves upstream with increasing Reynolds number. ThiB is not 
too surprising since the same behavior occurs with the upstream influ- 
ence. As the upstream influence increases, the "self-induced" pressure 
gradient will begin farther upstream and hence separation will occur 
farther upstream. For fixed Reynolds number, the upstream influence 
increased with decreasing Mach number, and accordingly (see fig. 9) 
the separation point, measured in multiples of S 0 *, moves Upstream with 
decreasing Mach number. 

Before closing the discussion on separation, the importance of D 
in the determination of the separation point should be discussed. The 
coefficients in the expression for A are Buch that, for fixed Reynolds 
number, Mach number, and deflection angle, the magnitude of A increases 
as D increases. The separation point has been seen to move upstream 
as A is increased in magnitude, hence an overestimation of the down- 
stream thickening would result in a slight overestimation of the upstream 
distance to the separation point. In figure 9 the variation of separa- 
tion point for .different Mach numbers has been plotted for fixed Reynolds 
number. The location of the separation point for the case D = 0 has 
also been plotted in this figure. These curves thus give the greatest 
lower bound of the separation distance since it is known that actually 
D >0. A comparison of results reveals that D has a very small effect 
upon the location of the separation point. For the higher Mach numbers 
in the range investigated, the percentage difference between the two cases 
may be fairly large, but the fact that the separation distance Is mea- 
sured in multiples of a boundary-layer displacement thickness must be 
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considered. In an actual experimental measurement of the absolute dis- 
tance to the separation pointy this difference will be extremely small. 
Hence , for practical purposes, one can consider D = 0 when estimating 
the separation point . The results obtained under such conditions will 
then yield a alight underestimation of the separation point. 


Pressure Distribution 


The pressure disturbance along the wall (fig. 10) decays exponen- 
tially from a def ini te value at the point of incidence to zero far 
upstream of that point. Downstream of the point of incidence the pres- 
sure rises to a maximum value before dropping to the value corresponding 
to reg ular reflection. This downstream overcompression has been observed 
experimentally, and it appears to be a characteristic feature of shock- 
wave interaction with a laminar boundary layer. The Lees' theory, as 
mentioned previously, failed to predict this downstream behavior. This 
is due to the fact that an incomplete solution for the pressure was used 
in the determination of the pressure distribution. The boundary layer 
was divided, longitud inally , into four regions, and it vhb assumed that 
the solutions to a third-order differential eguation were valid in each 
region. In the two regions that extended to positive and negative infin- 
ity, certain solutions could be rejected since they became infinite at 
the ends of their respective regions. In the finite regions, however, 
the complete solutions must be retained. The incompleteness of the Lees' 
theory, then, arises from the fact that only one term of the general 
solution was used in each of the finite regions. 

The linear theory yields pressure distributions that are s i milar 
irrespective 'of the size of the deflection angle. Now for very small 
angles, the experimental results exhibit the general behavior predicted 
by theory. In figure 10 the experimental values of the pressure distri- 
bution have been plotted for e = 1° and 1^, =2.05- (in this case, 

separation has probably not occurred since the wave Is quite weak.) The 
values of the pressure ratio were taken from figure 13 of reference 4 
nnri converted to the scale indicated. The upstream portion of this 
curve can be well represented by an exponential curve, thus verifying, 
at least for this case, the predicted exponential pressure rise. For 
larger deflection angles, the experimental pressure distributions are . 
characterized by the familiar pressure "bump" ahead of the point of inci- 
dence, thus ind icating that separation has occurred. The linear theory 
is inadequate, as regards predicting this upstream behavior; hence for 
large deflection angles the theory must be modified, possibly by taking 
account of the second -order effects upstream of the point of incidence. 
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SHOCK-WAVE INTERACTION WITH A TURBULENT BOUNDARY LAYER 
Method of Approximations 


When the flow in the boundary layer is turbulent, the Karman 
momentum- integral equation will take on the same form as for the laminar 
boundary layer except that the quantities involved are not the exact, but 
average, values. However, since general relationships between the shear 
stress and the mean velocity in turbulent layers have not yet been estab- 
lished, rigorous treatment of the problem at this time is, of course, not 
possible. Here, for the purpose of exhibiting the characteristic dif- 
ference between the laminar and turbulent cases, certain approximations 
are proposed. 

The momentum-integral equation expresses the shear stress at the 
wall in terms of the growth of the boundary-layer momentum thickness and 
the velocity gradient. Using experimental results as a guide, the rela- 
tive importance of these terms can be estimated. Experimental results 
of Fage and Sargent (reference 18) show that the shear stress in front 
of the shock is practically constant; behind the shock, it is very small. 
Therefore, unlike the laminar case, the shear stress to the first appproxi- 
mation, both in front of and behind the shock, can be regarded as con- 
stant and hence has no effect on the perturbed flow. It follows then 
that, in the case of a turbulent boundary layer, the growth of the momen- 
tum thickness is influenced primarily by pressure changes due to the 
presence of the, shock. 

To simplify the problem further, it is noted, for practical pur- 
poses, that the shape parameter H = 5*/0 is relatively insensitive to 
change even though there may be a considerable adverse pressure gradient. 

In the transonic case, where there is a normal shock in the local super- 
sonic region and a large change in H is anticipated, H at its maximum 
is only increased by a factor of about 1.2 (see reference 2). The reason 
for this is possibly the fact that the increase by the shock may be 
counterbalanced by a decrease due to compressibility effect. Moreover, 
in the momentum- integral equation, the coefficient of the velocity gra- 
dient is positive and usually greater than unity within the Mach number 
range; hence one can ignore the slight variation in H and consider it 
a constant. By the same token, the Mach number in this coefficient can 
also be considered to be constant. (The Mach number enters when the 
density is eliminated in terms of the velocity.) Consequently, the 
momentum- integral equation reduces to: 


d A 
dx 


a 3 


du* 

dx 


= 0 



bk 
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where 


a 3 = 2 + H - M^ 2 


Solution, for Outer Flow 

Similar to the previous case, for the upstream region 
= f ( 5 - m^T]). By applying the original condition dA/d| = 80/dq 

at t] = 0, there is obtained from the linear theory 


f . foe »W“3<S-"WI> + constant 

Since n^/a^ > 0 , the disturbances vanish far upstream. When cc.3 is 
evaluated the exponent mcc/ot^ is found to be about one order larger 
than A2 of the laminar case. This is the well-known experimental 

result that there is very little upstream influence in the turbulent 
case. This, therefore, confirms the hypotheses made in the previous 
section. 

For the downstream solution 


03 = sU - nwi) - + nyi) > 

The original differential equation remains the same; hence the down- 
stream solution would be 

„ BWWU-mcoT]) 26 . - ^ 

y'K ~ So e 3 £ + Constant 

J “oo 


It shows that if the solution 0i is continued to the downstream side 

the velocities would become infinite at positive infinity and must be 
rejected. There is therefore a principal difference between the laminar 
and turbulent cases. In the laminar case several solutions for the outer 
flow were obtained and the solution appropriate for either upstream or 
downstream could be chosen. In the turbulent case there is no choice 
since there is only one solution for the outer flow, which, if continued, 
fails at positive infinity. This Indicates that linearization of the 
flow is incapable of accounting for the flow near the point of incidence 
where the nonlinear effects become important. Since large changes in 
flow velocities can be brought about only through a shock., in the pres- 
ent problem a reflected shock must be considered. 
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From general considerations it can be argued that if one puts a 
shock wave between regions 2 and 3 and considers perturbations of a regu- 
lar reflection it may be possible to obtain a solution which will satisfy 
the conditions at infinity. By perturbing the regular reflection, there 
would be, in each region, undisturbed quantities plus their perturba- 
tions. Consider the total potential in region 3: 

$3 = |+ f(| - m 3 y) + g^| + n^y) 

At infinity f' = g* =0; along the reflected wave | = 0, to first 

order, the velocities are constant; hence 

f ± g* = Constant 

But to first order, along the wave, f'(0) = Constant; hence 

g'(S) = Constant. Therefore, g'(£) is constant throughout region 3» 

But g* = 0 at infinity; hence g* = 0 in region 3» 

Then at the boundary layer d A/d£ = -mof'd) where f is nondi- 

mensional, and the substitution into the momentum- integral equation 
yields 


-m 3 f' + ct^f’ ‘ = 0 

mo/ccof 

f = Constant + f Q e ->l J 

f' - f 0 | eV*' 

Now f 1 = 0 at infinity; hence one must choose f Q = 0, and thus 


f = Constant 


This means that behind the reflected shock, the flow, to the first order, 
is uniform. Since the flow behind the incident shock has undergone com- 
pression up to the reflected wave, it can be seen that, to satisfy the 
condition at infinity, the reflected shock must be followed by a very 
rapid expansion. Otherwise, the pressure after the second shock would 
have been higher than that after the regular reflection. It is there- 
fore concluded that in the case of a turbulent boundary layer the inci- 
dent shock is reflected as a shock and behind the reflected shock there 
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must be an abrupt expansion so as to cancel the overcompression brought 
about by the train of upstream compression waves. The existence of the 
sharp expansion behind the reflected shock was also confirmed by experi- 
ments (reference 9) • The pressure distribution and boundary- layer 
growth, to the first order, appear to be discontinuous as shown in fig- 
ure 12. 


Downstream Thickening for. Turbulent Boundary Layers 
The downstream thickening can be estimated by use of the formula 

r. v , ,2+H 

S3* Pl/ U l\ 

5 1* ” P3U13J 


Wow the difference between the laminar and turbulent cases is that H 
differs for the two cases. Since, in general, H-(- < H 1 , the downstream 

thicken in g in the turbulent case is less than the thickening in the 
corresponding laminar case. Empirical relations must be relied on to 
estimate H because of the lack of knowledge of the compressible turbu- 
lent boundary layer. According to Nitzberg and Crandall (reference 19), 
for local Mach numbers greater than the free-stream Mach number, but 
less than 1.4, the compressibility effect is well approximated by 

H = %0/i + § M 2 ) 


Use this' relationship to estimate the downstream thickening and 
then compare the results with the existing experimental data. By 
assuming a l/7-power law for the incompressible profile, 

H = 1.29^1 + O.^M^ 

The predicted thickening is shown in figure 13 . It shows that the down- 
stream thickening increases as the Mach number is decreased. This 
behavior is also present in the laminar case. 

Using this value of H, it was found that the predicted thickening 
for M ra = 1.44 and e = 4.5° is larger than the visually estimated 
thickening in the experimental case. On the other hand, the predicted 
thickening for Me = 2 and e = 6° is slightly less than the visually 
estimated thickening in the experimental case. It should be noted, how- 
ever, that, for M > 1.4, a relationship has been used that is supposedly 
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valid for M < 1.4. Moreover, this relationship merely accounts for the 
effects of compressibility. It would appear that some effect of the 
shock must also be included in the determination of H. However, at 
this time, owing to the lack of exact information, it must be neglected. 

Consequently, it appears that if the relation H = 1.29^1 + is 

used for Mach numbers not too much larger .than 1.4, one can expect to 
have a minimum estimate of the downstream thickening. 


SUMMARY OF RESULTS 


An investigation of the reflection of a weak shock wave from a 
boundary layer along a flat plate yielded the following results. 


Laminar Case 

>> 

1. In all cases investigated, the pressure along the wall overcom- 
presses downstream of the point of incidence. The pressure disturbance 
decays exponentially from a definite value at the point of incidence to 
zero far upstream of the point of incidence. Downstream, the pressure 
rises to a maximum value and then falls, gradually, to the constant 
value corresponding to regular reflection. 

The exponential pressure rise appears to be verified in the case of 
a shock deflection angle of 1°, since separation has probably not 
occurred. For larger deflection angles (the next larger angle for which 
there are experimental pressure distributions is 3°), the experimental 
pressure distributions exhibit the familiar pressure bump between the 
separation point and the origin. For these angles, a true comparison 
between experimental results and theoretical results cannot be made 
since the present theory does not account for the effects of separation. 

2. If the upstream influence is considered to be the distance to 
the point at which the disturbance has decayed to a specified fraction 
of its amplitude at the origin, the upstream influence, when measured in 
multiples of the boundary-layer displacement thickness 5 0 *, is found to 

increase with increasing Reynolds number: x/5 0 * « ^Re/m^^/ 2 where x 

is the coordinate para! lei to the flow direction. Re is the Reynolds 

number, = \ju oo 2 - 1, and Moo is the free-stream Mach number. For 

decreasing Mach number, the upstream influence also increases. If the 
disturbance is considered to decay to, say, 5 percent of its amplitude 
at the origin, the upstream influence for Mo, ~ 2 and Re ~ 1500 is of 
the order of 30 boundary-layer displacement thicknesses. 



NACA TN 2868 


kQ 


3- The pressure gradient is such that the boundary layer may separate 
ahead of the point of incidence. In the event of separation, an increase 
of shock strength, for fixed Mach numbers and Reynolds number, increases 
the distance between the separation point and the point of incidence. 

For fixed Mach number and shock-deflection angle, the separation point ' 
measured in multiples of 5 0 * moves upstream with increasing Reynolds 
number. For the case of M^, =.1.^(4 and Re = 2000, the boundary layer 
separates for all flow-deflection angles e > 1 . 03 °. 

For fixed Reynolds number and shock-deflection angle, the separa- 
tion point, measured in multiples of 5 0 *, moves upstream with decreasing 
Mach number. 

For a complete determination of the constants of integration, an 
estimate of the downstream boundary-layer thickness was required. For 
the cases investigated, the effect of the downstream thickening on the 
outer flow is rather small. In fact, for practical purposes, this effect 
of downstream thickening can be neglected in the determination of the 
separation point. The distance between the origin and the separation 
point will then be slightly underestimated. 

4. The present theory is applicable only in the case when a weak 
shock is incident upon a laminar boundary layer. In addition to the 
outer flow field, the boundary- layer displacement thickness has also 
been linearized. The latter linearization enables a linear differential 
equation with constant coefficients to be obtained for the perturbation 
velocity potential. This equation immediately yields the general solu- 
tion of the outer flow. Upon investigating the size of the perturba- 
tions, the maximum velocity perturbations are found to be about 10 per- 
cent or less. The thickness perturbations, on the other hand, are much 
larger, being about 30 percent near the origin; thus the linearization 
of the displacement thickness becomes questionable as the point of inci- 
dence is approached. Had the displacement thickness not been linearized, 
a rather complicated nonlinear differential equation would have been 
obtained from the boundary condition for the perturbation velocity poten- 
tial. A solution of this equation would be expected to yield more accu- 
rate results. However, it is problematical as to whether the equation 
could be solved without the imposition of additional assunqrtlons which, 
in themselves, might nullify- any accuracy that the nonlinear boundary 
condition may provide. 


Turbulent Case 

In the case of shock-wave interaction with a turbulent boundary 
layer the upstream influence is found to be considerably less than in 
the laminar case. In addition, to first order, the incident wave must 
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"be reflected as a compression wave followed immediately by an expansion 
wave, so that the end pressure condition is satisfied. 


Cornell University 

Ithaca , N. Y., January 11, 1952 
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Figure 3»- Displacement -thickness ratio 
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(a) Various values of M^. 

Figure 1)-. - Measure of upstream influence. logg (l/a)j 

disturbance at x<j/ S 0 * is a percent of disturbance at origin. 
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Figure 6.- Reynolds number against Re = U5 q*/v 0 . 





Figure 8.- Separation point against deflection angle. = 1 . 44 j 

Re « U5 o */v 0 . 
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Figure 11.- Perturbed "boundary-layer displacement thickness. = 1.44j 

Re = U8 0 */ v o = 500- 
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